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I. INTRODUCTION 



The supersymmetry (SUSY) algebra in quantum mechanics (QM) began with the work 
of Nicolai [1] and was elegantly formulated by Witten [2]. This approach has attracted 
interest and was applied to construct the spectral resolution of solvable potentials in many 
areas of physics [3]. SUSY QM was also formulated by Gendenshtein who used the shape 
invariance property [4] . The hydrogen atom was studied via SUSY QM in the non-relativistic 
context by Kostelecky and Nieto [5]. They used the SUSY QM for spectral resolution and 
also for calculating transition probabihties for alkali-metal atoms. Also, Zhang et al. have 
considered interesting applications of a semi- unitary formulation of SUSY QM [6] . 

The Dirac-Coulomb problem is an exactly solvable problem in relativistic quantum me- 
chanics and the solution can be found in books on quantum mechanics, see for instance [7]. 
The Dirac-Coulomb problem has also been studied via SUSY QM [9-17]. Our purpose in 
this paper is to obtain the complete energy spectrum and the energy eigenfunctions of the 
Dirac-Coulomb problem using the shape invariance and new generalized ladder operators. 
The Lie algebra associated with these operators for the shape invariant potentials has been 
presented by Fukui-Aizawa [18] and Balantekin [19]. This formalism has been applied to 
the exactly solvable potentials in non-relativistic quantum mechanics [20,21]. 

It is particularly simple to apply SUSY QM to shape-invariant potentials because their 
SUSY partners are similar in shape and differ only in the parameters that appear in them. 
More specifically, if V-{x;ai) is any potential, adjusted to have zero ground state energy 
= 0, its SUSY partner V+i ) must satisfy the requirement V+{x;ao) = V^{x;ai) + 
R{ai), where is a set of parameters, ai a function of the parameter oq and R{ai) is 
independent of x. In this case, one can determine the energy levels for V+{x;ao), which 
correspond to = J2i=i R{(ii)- In this context, recently, some relativistic shape invariant 
potentials have been investigated [22]. 

The SUSY hierarchical prescription was used by Sukumar to solve the energy spectrum of 
the Dirac-Coulomb problem [10]. In this work the Fukui-Aizawa-Balantekin [18,19] approach 
to the Dirac-Coulomb problem is investigated using SUSY QM. 

The paper is organized as follows. In section II we realize a graded Lie algebra struc- 
ture in terms of the 4x4 matrix supercharges, analogous to Witten's SUSY algebra for the 
Dirac radial equation associated with the hydrogen atom. In section III, using the shape 
invariance condition wc deduce new generalized ladder operators in relativistic quantum me- 
chanics, via supersymmetry, in order to build up the energy eigenvalue and eigenfunctions 
of super symmetric partner potentials. In section IV some concluding remarks are given. 



II. THE DIRAC-COULOMB PROBLEM AND SUSY 



In this section, we adopt the Sukumar's approach to construct the 4x4 matrix super- 
charges. The Dirac radial equation for the hydrogen atom can be written as 

0^ , ^ ( k \ ( G \ _ ( Q ai\(G\ 




(1) 



where k is an eigenvalue of the Dirac operator K = /3(E ■ L -f 1), 7 = ai — m + 

E, a2^m-E, \ k \^ j + ^ (A; = ±1, ±2, ±3, • • •), E = (Ei, E2, S3) are a set of the 
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Pauli matrices, L is the angular momentum operator and 1 is the 2x2 unity matrix. Now, 
let D be the operator which diagonalizes the matrix that appears in the interaction term 

D-\kaz-'na2)D = sa^, (2) 

where s = (A;^ — 7^)^. Then we get 

(s + A;)l + 7(7i. (3) 



Thus, we obtain 



where 



The eigenvalue equations ior k —\k \ and k — — \ k \^ respectively, become 

(\ k I m\ ~ f d s j\ ~ ,„ , 



s El \dp p s 



' - - G+ = -- + - - 4 F+, (6b) 



s y \ dp p s 



^-B\F-=[-r + ---\G-^ (6c) 



s y \dp p s 



= + (6d) 



s £^ y y dp p s 



where F± = F(± | k |) and G± = G{± \ k |). 
Defining the intertwining operators as 



^0 
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where I denotes the 2x2 unit matrix and 

m 

E"' ' " s 



O = -(Ti + — '-a2, (9) 



we get 



4'Hf)=o(f], (10) 



where I denotes the 2x2 unit matrix. 
Using the result 

Ao = -(7,A+a,, (12) 
we see that there exist the supersymmetric partner eigenvalue equations: 

and 

The mutually adjoint non-Hermitian supercharge operators for Witten's model are given 

by 



so that the SUSY Hamiltonian H takes the form 



// = [g+,Q_]+= ( ^- ^ ), (16) 



and satisfies [i?, Q±\- — 0. 

The pair of SUSY hamiltonians is given by 



= + (17) 



dp^ dp 

where the matrix superpotential is given by W{p) — ~ s) c^a- 
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III. SPECTRAL RESOLUTION VIA LADDER OPERATORS 



Let us now build up the energy eigenvalues and eigenfunctions of supersymmetric partner 
potentials using the shape invariance condition and the generalized ladder operators. Prom 
the last section one obtains the following matrix forms of the pair of SUSY potentials: 

/ sjs-l) _ 27 , 7^ Q \ 

V-{p.\s)^[ ^ , ^ , (18) 

^+(P,A,5)= ^ £M_27 , 7^ • (19) 

Although the SUSY partner potentials V(±) are not shape invariant, we can see that their 
components are: 

V'(+)ii(p, 7, s) = V'(-)ii(p, 7, s + 1) - ^gj^Yf ^ ^' ^^^^ 

2 2 

y(_)22(p,7,*) = 1^(+)22(P,7,S + 1) - + ^- (21) 



From (26), (20) and (21) one can written 



72 ^2 ^2 ^2 



ii„(a,) = fl.W = -^ + ^ = ^-^. (22) 



v2 ^2 ^2 



where = s and ai = Oq + 1, so that = ao + i, -Rii(aj) = , 2-\i ^ \ — j-^ + Thus 
we get the following energy eigenvalues of = -f^(+)22 : 

n ™ / 1 1 \ _'v2 n/2 

^tV(«0' «l-iy (s + n)2 s2 

Comparing (23) with the eigenvalue equation (14) we obtain the energy eigenvalues of the 
hydrogen relativistic atom, viz., 

7^ J7j2 ^2 ^2 

1 + ~ E^Cnia = ~2 ~ / „ , „\2 ' (24) 



S 



^;(n)2 52 (5 + ^)2 ' 



providing 



£;(n) 



m2 



\ 1 + ^ 



-, ri = 0,l,2,---, (25) 



which is in agreement with the result obtained by Sukumar using the SUSY Hamiltonian 
hierarchy method [10]. 



5 



Note that the shape invariance condition is associated with translation of the parameters 
a's, so that the Eq. (20) can be written in the following form 

4i Vo)4tVo) = A[t\a,)A[-\a,) + i?n(ai), (26) 

where 

4f(s)-±:^ + ---, ao = s. (27) 
dp p s 

Following Fukui-Aizawa-Balantekin approach [18,19], we obtain the following ladder op- 
erators 

B4s)=T\s)A^,{s), B+{s) = bUs), (28) 

where T(s) being a translation operator defined by 

T{s)=e^, (29) 

with T^{s) = e~^. In this case we have the identity i2(a„) = T(ao)-R(a„_i)T^(ao) and the 
folowing algebra 

[B-, B+] = T{ao)R{ao)T^{ao) = i?(ai), i?(a,)S+(ao) = S+(ao)/?(a„_i), (30) 

where the shape invariance provides us the translations of the parameters a„, viz., a„ = 
Gq + n, valid for any n. Thus, it is easy to see that the operators i?±(ao) and R{an) satisfy 
the following commutation relations 

[i/(_)n, Bl] = (i?(ai) + R{a2) + ■■■ + i?(a„))S;, (31) 

and 

[i/(_)n, B^] = -Bl{R{a^) + Ria^) + ■■■ + R{a^)), n = 1, 2, • • • , (32) 

with = B^B_ = A^^^A['^\ Consequently we see that the component eigenfunc- 

tion of the n-th excited stated, is given by 

f!") oc S!^(s)gL°^(/9;s), n = l,2,3,---. (33) 

The ground state eigenfunction must be annihilated by -B_(s), then 

A-is)G^^\p;s) = Q, (34) 

which leads us to the following physically acceptable solution 

GL°^(p;s) = iVGA->, (35) 

where Nq being the normalization constant. 

^From (28) and (29) we see that the raising operator may be written as 
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d S j\ a 
— + - - - es^. 
dp p s J 



(36) 



Consequently, for the first excited state one may write 



F^'\p;s)=B^{s)G^^\p;s) 




s+l — ^ 



(37) 



Finally we would like to call attention to the fact that the above formalism may be ap- 
plied to the exactly solvable potentials in relativistic quantum mechanics [22,23]. Also, the 
applications of semi-unitary transformations to construct supersymmetric partner Hamilto- 
nian in non-relativistic quantum mechanics [6] can be implemented for the Dirac-Coulomb 
problem. 



In this paper we investigated the Dirac-Coulomb problem via supersymmetry in quantum 
mechanics. The shape invariant formalism for the supersymmetric partners is applied to 
obtain the complete energy spectrum and eigenfunctions of the Dirac-Coulomb problem. 

Our approach uses the algebraic structure for shape invariant potential, recently proposed 
by Fukui-Aizawa-Balantekin [18,19] . This approach is different from the SUSY Hamiltonian 
hierarchy method applied by Sukumar [10]. 
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